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SUMMARY , 

h-f?- 6 7 

The various hypersonic flow parameters for helium have been estimated by use 
of the Beattie-Bridgeman equation of state for helium at stagnation pressures up 
to 50,000 lb/sq in. abs and stagnation temperatures up to 10,000° F. Calcula- 
tions were limited to densities below the critical density. The results were 
obtained to establish reasonable estimates of the departure from ideal-gas behav- 
ior over a broad practical range of stagnation pressure and temperature at Mach 
numbers up to 100. The results are presented as correction factors which are to 
be applied to ideal-gas parameters for a specific-heat ratio of 5/3* 


INTRODUCTION 


Hypersonic helium tunnels have been used to study fluid dynamic problems for 
a number of years. Most of this work has been carried out at Mach numbers in the 
neighborhood of 20 with stagnation temperatures not much different from room tem- 
perature and stagnation pressures up to about 3>000 lb/sq in. abs. For these 
stagnation conditions, there can be appreciable errors in the determination of the 
various flow parameters when based on the ideal-gas parameters for a specific-heat 
ratio of 5/3 because of deviations from ideal- gas behavior at the higher densities 
For example, results of calculations presented in reference 1, which are based on 
the thermodynamic properties of helium reported in reference 2, indicate that the 
free-stream temperature, pressure, and density are approximately 5> 9> and 15 per- 
cent greater, respectively, than the ideal-gas values for a specific-heat ratio of 
5/3 when the stagnation temperature is 80° F and the stagnation pressure is 
3,000 lb/sq in. abs. At the same temperature but higher stagnation pressures, 
the errors are even greater. 

Reference 1 presents the results of calculations as correction factors which 
must be applied to the tabulated ideal-gas parameters for a specific-heat ratio 
of 5/3* The ideal-gas parameters for helium which must be used in conjunction 
with these correction terms are presented in reference 3 for Mach numbers up to 
4-0 and in reference 4- for Mach numbers from 4-0 to 100. 

The thermodynamic data for helium presented in reference 2, which were used 
in calculating the correction factors in reference 1, are limited to the pressure 



range of lL .7 to 6,000 lb/sq in. abs and temperatures up to 600° F. Likewise, 
the correction terms of reference 1 are limited to stagnation pressures up to 
6,000 lb/sq in. abs and temperatures up to 600° F. 

It appears possible to operate a hotshot-type tunnel, with an appropriate 
arc chamber, with helium as the test gas at Mach numbers of 60 or so at stagna- 
tion pressures as high as 50,000 lb/sq in. abs over a broad range of stagnation 
temperatures. In addition to this approach to obtaining higher Mach number 
flows in helium, work at the Langley Research Center is in progress to develop 
a conventional Mach 40 helium blowdown tunnel capable of operating at stagnation 
pressures up to 15,000 lb/sq in. abs with a stagnation temperature of about 
600° F. 

Because of the limited pressure and temperature ranges of the available 
thermodynamic properties of helium, the correction terms presented in reference 1 
fall short of the stagnation conditions for operating the previously mentioned 
facilities. It is the purpose of this present work to obtain reasonable esti- 
mates of the correction terms for the hypersonic flow parameters for helium at 
stagnation temperatures up to 10,000° F and pressures up to 50,000 lb/sq in. abs 
at Mach numbers up to 100 or that limited by equilibrium condensation of helium. 


SYMBOLS 


a velocity of sound 

a',A 0 ,b,B 0 ,C constants in equation (l) 


Cp specific heat at constant pressure 

c y specific heat at constant volume 

c v specific heat at constant volume for ideal-gas conditions 

F correction factor (The correction factor is a particular flow parameter 

calculated by dividing the real-gas value by the ideal-gas value for 
7 = 5/3 at the same Mj_. For example, F pjl = pq/p-^q where both 

Pq and Pq^q are evaluated at the same Mq. ) 

h specific enthalpy 

M Mach number 


p pressure 

q dynamic pressure 

R gas constant 


2 


s specific entropy 

T temperature 

V velocity 

W limiting velocity ratio 

Z compressibility factor 

7 ratio of specific heats 

p density 

Subscripts: 
c 

1 

t,l 
t,2 
1 

2 

THERMODYNAMIC RELATIONS 


conditions for equilibrium condensation 
based on ideal gas and 7=5/5 
stagnation conditions upstream of shock 
stagnation conditions downstream of shock 
free-stream conditions upstream of shock 
free- stream conditions just downstream of shock 


Equation of State 

The thermodynamic properties of helium presented in reference 2 include 
pressures up to 6, 000 lb/sq in. abs and temperatures up to 600° F. For the most 
part, these data are based on the Beattie-Bridgeman equation of state. In order 
to calculate the desired flow parameters for stagnation conditions beyond this 
limited range, the thermodynamic properties of helium must be extended. In the 
present work the Beattie-Bridgeman equation with appropriate constants for helium 
was assumed to be an adequate representation of the equation of state for helium. 

It should be emphasized that results obtained by the employment of this equation 
are to be treated only as first-order estimates since the constants of the equa- 
tion are determined for a rather limited range. This equation is not reliable 
above the critical density so it was necessary to restrict the use of it to den- 
sities below the critical density. The critical density of helium is lb/cu ft. 
The Beattie-Bridgeman equation is 

p = RT(1 - e)(l + Bp) p - Ap 2 (1) 

where A = A^l - a'p)j B = B Q (l - bp)j e = Cp/T^. The value of the constant b 
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for helium is zero so that equation (l) can be rewritten in terms of the con- 
stants a', Aq, Bo, and. C. Thus, 

p = AoP 2 (a'p - 1) + Rp(BoP + 1)T - RCp 2 (BoP + l)-^ ( 2 ) 

ir 

This equation is used, in conjunction with various thermodynamic relations to find 
the expressions for enthalpy, entropy, and velocity of sound, each as a function 
of density and temperature. 


Enthalpy 

The expression for enthalpy is found by beginning with the relation, 


dh 



dT + 



(3) 


From the definition of enthalpy and c v , the first term on the right-hand side of 
this equation can be written as 



00 


The second term of equation (3) can be written as a function of T, p, and p 
by introducing the equation. 


dh = T ds + — dp 
P 


which can be rewritten as 



and by using the Maxwell relation, 



Thus, 



-T/dp\ 

P 2 W p 


+ -(^) 
P\dp/T 


(5) 


( 6 ) 


(7) 


( 8 ) 


Since h is a point function, equation (3) may be integrated over a two-step 
path by first holding p constant and integrating with respect to T and then 




by holding T constant and integrating with respect to p. Before this integra- 
tion is performed, p is eliminated from equations (4) and (8) by use of equa- 
tion (2). By direct differentiation of equation ( 2 ), 

(tf)p ' M** + x ) + + (9) 

and 

(i) T ' 1,1 + 2 ( rb3t - 1 - *°y + 3 ( A ° a ' - wy <io) 

These two equations are then used in equations (4) and (8). The density level 
chosen for the integration along the constant density path from a reference tem- 
perature to a desired temperature can be made as small as desired so that c v in 
equation (4) can be replaced by c v *, which is the value of c v as p -» 0 and is 
equal to a constant of (3/2)R. The second integration is carried out along the 
path of constant desired temperature from this very low reference density up to 
the desired density. Since the reference value of the enthalpy is an arbitrary 
constant, it can be set equal to zero at the reference conditions. Upon integra- 
tion, the expression for enthalpy as a function of p and T becomes, 

h = (c v * + r)t - ^2Ao - RBqT + 4 pV - - 3Aoa'jp 2 (ll) 


Entropy 

The expression for entropy as a function of T and p is found by beginning 
with the relation. 



The first term on the right-hand side of this equation can be written as 



(13) 


and the second term is given by equation (j). Upon substitution of equations (7) 
and (13) into equation (12) and integrating over the same two paths for which the 
enthalpy was determined and setting the value of entropy equal to zero at the 
reference conditions, the expression for entropy as a function of T and p is 


s = c v log e T - R 


log e p + ^Bq + pjp + p : 


( 1*0 
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Velocity of Sound 


The expression for the velocity of sound is found by beginning with the 
relations, 


a 


2 




( 15 ) 


and 


C P 





(16) 





(17) 


Equation (l6) is divided by c v and multiplied by the quantity 

to equation (17). Then to this resulting equation, the quantity 
to give the desired equivalent of equation (15), 



according 
is added 


2 _ /dp\ T lb p\ 2 

WT + p 2 c v Wp 


(18) 


The two quantities, 


'bp\ . 




P/T 


, are given by equations (9) and (10), 


respectively. The expression for c y at a given temperature can be determined 
at a desired density p from the equation. 


= c v * + 


r (— i 

J o \8p /t 


dp 


(19) 


where c v * is the heat capacity at constant volume for the condition of p -» 0. 
The integrand of equation (19) is replaced by the relation. 



T f 6 2 p\ 

p 2 Ut 2 / p 


and by differentiation of 


equation 


(9), 



p 


- - 6bc p 2 (^p + 1 )^. 


( 20 ) 


( 21 ) 
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The integration of equation (19) with the substitution of equation (2l) yields 

c v = c v * + 3 p(BqP + 2) (22) 

which is used to evaluate c v in equation (l8). As shown in the next section, 
these thermodynamic relations are used to determine the various flow parameters 
for helium for a number of stagnation conditions and Mach numbers. 


CALCULATION OF FLOW PARAMETERS 


The foregoing derived relations for enthalpy, entropy, and velocity of sound, 
along with the equation of state, are used to calculate the various flow param- 
eters. The numerical results for a range of stagnation conditions and Mach num- 
bers were obtained with an IBM 7090 electronic data processing system. In this 
section the method of computation is discussed. 

The first step of the calculation scheme is to establish the stagnation con- 
ditions by specifying the stagnation temperature T^ ^ and stagnation density 

The stagnation pressure p^ stagnation enthalpy h^. and stagnation 
entropy s t,l are "then calculated by application of equations (2), (ll), and 
(l4), respectively. 


Free- Stream Conditions 

Various free- stream conditions are determined for an arbitrary set of stag- 
nation conditions by choosing a series of free- stream temperatures T^_, which 

are successively less than For each value of Tj_, the value of the free- 

stream density p-^ is determined for an isentropic expansion from the stagnation 
condition by setting s in equation (l4) equal to and T equal to T-^ and 

iterating for p^. After p^ is determined for a given value of T-^, equa- 
tions (2) and (ll) can be used to find p^ and hj_. The corresponding velocity 

of sound is determined from equation (l8) with use of equations ( 9 ), (10), and 
(22). The free-stream velocity is then determined from the energy equation, 

Vl = )j 2(h t ,i - hx) (25) 

and the free-stream Mach number M-j_ is simply found by dividing Vj_ by a-j_. 

To this point, for a specific stagnation condition, Tq_, p^, p]_, Vj_, and M^ 
are determined. A number of values for T^ can be taken to find these same 
parameters for a range of M^_. 
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Conditions Just Downstream of Normal Shock 


The conservation relations for mass, momentum, and energy are 


Pl V l = 


p 2 v 2 


( 24 ) 


P]_ "*■ PqVp — P2 P2^2 


(25) 


and 


\,1 - + | V l 2 - *2 + | V 2 2 < 26 > 

and are used in conjunction with equations (2) and (ll) to determine the condi- 
tions just downstream of the normal shock. This determination is accomplished by 
eliminating V 2 from equations (24) to (2 6) and replacing p 2 by equation (2) 

and hp by equation (ll) followed by a double iteration for p 2 and T 2 . Once 
p 2 and T 2 are determined, p 2 is calculated from equation (2), from 

equation (ll), s 2 from equation (l4), V 2 from equation (23), and from 

equation (l8). 


Stagnation Conditions Downstream of Normal Shock 

The stagnation conditions downstream of a normal shock must satisfy two con- 
ditions, Sp 2 - s 2 and - &tj2 = kt,l* Since s 2 and h-^q have been calculated 

for a given stagnation condition and free- stream condition, &t,2 an< ^ ^t,2 are 

known. The simultaneous solution of equations (ll) and (l4) by double iteration 
determined p^ 2 and T-^ 2 . The value of P-t^2 is then determined from equa- 
tion (2). 

The various flow parameters are therefore determined for a number of combi- 
nations of stagnation conditions and Mach numbers based on the Beattie-Bridgeman 
equation for helium. The constants used in this work which appear in the Beattie- 
Bridgeman equation, according to reference 2, are 

a' = 0.23963 




5.0906 


ft3\ 2 ^f_ 

W in. 2 
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B q = 0.056063 

■ LD m 


' 93lKl7 s ~° r3 


R 


2.6829 


( ft V lb m)( lb f/ ln - g ) 

°R 


where the subscripts m and f indicate mass and force, respectively. 

In these calculations the various parameters were made nondimens ional by 
dividing by an appropriate reference condition. The following parameters were 
obtained for a number of stagnation conditions and various values of M^: 

T l/ T t,l> Pl/Pt,l> Pl/Pt,l> 9 . 1 /Pt, 1’ w l> t 2/ t 1> P2/Pl> P 2 / Pl> m 2> 

T t 2 /^t 1> p t 2 /Pt \> and p t 2 /Pt 1' For ease of P resentation > these param- 
eters were divided by the corresponding ideal-gas parameter with y = 5/5 for 
the same value of Mi to give the correction factors. For example, the cor- 
rection factor for the free-stream temperature, j, is ( T i/ T t, l)/( T l/ T t, l) 

or for a given value of M-j_. 


Vapor- Pres sure Curve for Helium 

In order to calculate the maximum Mach number which can be obtained from an 
equilibrium isentropic expansion in helium without helium condensation, the 
vapor-pressure relation presented in reference 5 is used. For T c > 3*9^2° R 

log 1Q P c = - + 2,208 1Og 10 T c " 0-06027 ( 27a) 


and for T c < 3>9 1 +2 0 R 


log 10 P c = - 6 ‘j^ 62 + 0.922 log 10 T c + 1.0860 


(27b) 


where T c is the temperature at which helium just condenses when the pressure for 
equilibrium conditions is p c . This relation is applied to find the Mach number 
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at which helium condensation just occurs for equilibrium flow for a given set of 
stagnation conditions by first calculating the free- stream pressure and Mach 

number Mq_ for a particular value of the free- stream temperature Tj_. This value 
of Tp is used in conjunction with equation (27a) or (27b) to find the pressure 
at which equilibrium condensation could occur. If this calculated pressure p c 
is greater than the free-stream pressure p^, condensation will not occur. On 
the other hand, if p c is less than pp, condensation will occur for equilibrium 
flow. The state such that p^ = P c at the same time T-^ = T c is the condition 

for maximum Mach number without condensation. The value of Mach number at which 
helium condensation will just occur for a given stagnation condition was deter- 
mined by an iterative procedure and was the maximum value used in these 
calculations. 


RESULTS AND DISCUSSION 


The results of the calculations for helium based on the Beattie-Bridgeman 
equation of state are presented in figures 1 to 8 as correction factors to be 
used in conjunction with the ideal-gas parameters for helium presented in refer- 
ences 3 and 4. A comparison between the values of the correction factors for 
hypersonic flow parameters in helium presented in reference 1 and those calcu- 
lated in the present work shows an agreement that is always better than 3 percent 
even when the effect of Mach number is considered. As was noted previously, the 
work of reference 1 was based on the assumption of ideal-gas behavior below 
1 atmosphere, which would cause the various correction factors to become inde- 
pendent of Mach number so long as the expansion was carried below 1 atmosphere. 

In the present work the Beattie-Bridgeman equation was applied over the entire 
flow process and, as a result, a small dependency on Mach number was noted even 
for expansions carried below 1 atmosphere. This small difference between the 
results of reference 1 and the present work is attributed to the slightly dif- 
ferent methods of calculation. 

The correction factors for Wp, Pp/Pi* ^ 2 . 1 ^ 1 ’ an<i ^2 are n °b presented 
since they were all within 1 percent of unity for the entire range of calculation. 

The results presented in reference 1 were limited to stagnation pressures up 
to 6,000 lb/sq in. abs and temperatures up to 600° F, whereas the present work is 
for stagnation pressures up to 50,000 lb/sq in. abs and temperatures up to 
10,000° F with the restriction that the calculations are limited to densities 
below the critical density. Because of the broad range of conditions considered, 
the results presented herein are to be considered simply as reasonable estimates 
of the correction factors for the various hypersonic flow parameters for helium. 

The use of the correction factors presented in figures 1 to 8 is easily shown 
by an example. Let us assume that a hypersonic helium tunnel operates with a 
stagnation temperature of 600° F and a stagnation pressure of 10,000 lb/sq in. abs 
at a ratio of pitot pressure to stagnation pressure of 3-50 x 10“^. For these 
stagnation conditions, the correction factor Fp ^ ^ i- s approximately 1.185 
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(see fig. 7) with some uncertainty due to Mach number. The ideal-gas value for 


- 1 is therefore ( 3-50 X 10-^/1.185, 


>-4 


The actual Mach num- 


or 2.95 X 10" 

her for these conditions is then found from reference 4 to be approximately 42. 60. 
Since there was initially a question as to which Mach number to use in reading 
figure 7, the value of should be read again with the added knowledge that 

Mq is very close to 42. 60. Now by considering the small effect of Mach number in 
figure 7, the value of 

Pt 2 /Pt 1 of 2.96 X 10-4, so that the Mach number is 42.55* 


F_ , o is l.l8l, which gives an ideal-gas value for 
p, x., c; 


Now that the value of Mq is determined, the various ideal- gas parameters 
can be found for Mq = 42.55 from reference 4 and the corresponding correction 
factors obtained from figures 1 to 8. For example, the actual value of T-^^T^. 

is determined by multiplying the ideal-gas value of I .65 X 10“3 found in refer- 
ence 4 for M^ = 42.55 by a correction factor of 1.102 found from figure 1 for a 

stagnation temperature of 600° F and a stagnation pressure of 10,000 lb/sq in. abs 

which gives 1.102 ( 1.65 X 10“3), or 1.82 X 10”3. The other flow parameters are 
determined in the same manner. 


Figure 9 shows the maximum Mach number that can be attained by the isentropic 
expansion from various combinations of stagnation temperatures and pressures for 
equilibrium conditions. The values presented in this figure are also to be con- 
sidered simply as reasonable estimates. 

Finally, figure 10 shows the compressibility factor Z plotted as a function 
of pressure for a number of temperatures as calculated from equation (2). For a 
given temperature and pressure, the density can be found from the relation 


p = ZpRT 


( 28 ) 


and figure 10. 


CONCLUDING REMARKS 


The various hypersonic flow parameters for helium have been estimated by use 
of the Beattie-Bridgeman equation of state for helium for stagnation- pressures up 
to 50,000 lb/sq in. abs and stagnation temperatures up to 10,000° F. Calculations 
were made for Mach numbers up to 100 or the Mach number corresponding to the point 
of equilibrium helium condensation, whichever was smaller. These calculations 
were restricted to densities below the critical density and are considered to be 
only reasonable estimates of the flow pa 'ameters for helium. The results are pre- 
sented as simple correction factors which can be used in conjunction with the tab- 
ulated ideal-gas parameters presented in NACA TN 4063 and NASA TN D-1252 for a 
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specific-heat ratio of 5/3- This work extends the stagnation pressure and temp- 
erature range considered in NASA TN D-462. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va. , January 21, 1963- 
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Correction facto r, 



50 x 10 3 


Stagnation pressure, p t p, lb/ sq in. abs 


Figure 1.- Correction factor F_ as a function of stagnation pressure for various stagnation 

i. f -L 

temperatures and Mach numbers . 
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Correction factor 










Correction factor . 



Figure 3. - Correction factor 


Fp^i as a function of stagnation pressure for various stagnation 
temperatures and. Mach numbers. 
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Stagnation pressure, p t lb/sqin. abs 


Figure 4.- Correction factor F q as a function of stagnation pressure for various stagnation 

temperatures and. Mach numbers . 





Correction factor . 





Hi 

m 

m 

M 

m 

m 


M 

w 

m 

m 

affi 

m 

BS 

«b 

m 



m 

m 

im 

a 

M 

oB 

a 

m 

■? 

n 

m 

a 

a 

a 

a 

m 

a 

a 

a 

all 

a 


a 

K 

a 


a 

a 



a 

t i 


a 




a 

■SB 

!l!!ffi!li8iiiiii 

j 

1 

a 

a 

i 

a 

a 

i 

i 

1 

a 

a 

i 

a 

a 

1 

1 

| 

i 

a 

a 

i 


i 

1 



1 




a 



1 






Figure 5-- Correction factor F T ^ 2 as 8 function of stagnation pressure for various stagnation 

temperatures and Mach numbers. 
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Figure 6.- Correction factor as a function of stagnation pressure for various stagnation 

temperatures and Mach numbers. 
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Correction factor, F p>tj 2 



Figure 7- - Correction factor 


Fp, t,2 as a function of stagnation pressure for various stagnation 
temperatures and Mach numbers. 
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Correction factor. 



Figure 8. - Correction factor 


Fp t 2 as a Function of stagnation pressure for various stagnation 
temperatures and Mach numbers. 
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Figure - Mach number at which condensation just occurs for equilibrium conditions as a function of stagnation pressure 

for various stagnation temperatures. 



Compressibility factor, 



22 


NASA-Langley, 1963 L— 3^6l 


